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METHOD AND DEVIC E FOR CORRECTING PROXIMITY EFFE CTS 

The present invention relates to a method and 
device for determining the exposure dose of an electron 
beam required per pattern position to obtain a desired 
pattern in a coating. 
5 In the manufacture of the latest generations of 

integrated circuits use is preferably made of focused 
electron beams in lithographic processes instead of 
making use of the usual optical lithographic techniques, 
since these latter techniques are subject to limitations 
10 in terms of resolution as a result of diffraction of the 
used laser light. The resolution of the , integrated cir- 
cuit obtained with such electron beam lithography is 
greater, although it is limited by scattering of the 
electrons in the coating. Methods are known for minimi - 
15 zing scattering effects or compensating therefor in 
advance and thereby improving the resolution of the 
-obtained integrated circuits. 

The known methods have the drawback however 
that scattering effects can themselves only be minimized 

2 0 to a limited degree, while advance compensation according 

to the known method requires many calculations and there- 
fore needs a long calculation time. For the manufacture 
of integrated circuits for instance a very large number 
of pattern points, often in .the order of magnitude of 10 
25 pattern points, must be •written", while the number of 
calculations required for this purpose amounts to a 
multiple thereof. As a result a practically real time 
precompensation for the smearing or blurring effects 
cannot be implemented. 

3 0 The object of the present invention is to 

obviate this drawback and also provide additional advan- 
tages . 

The present invention therefore, relates to a 
method for determining the precompensated pattern of 
3 5 exposure doses of an electron beam required per pattern 
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position to obtain a desired pattern in a coating on a 
substrate, comprising of: 

- determining the smearing function (blur 
function) of the electron beam; 
5 - determining the precompensated pattern with 

the smearing function and the desired pattern, wherein 
the determination is performed such that the exposure 
doses contain almost exclusively positive values and that 
the exposure doses are at least to some extent smooth 
10 relative to each other. 

Since a negative value for the exposure doses 
of an electron beam has no physical significance and 
cannot therefore be realized, the determination of the 
exposure doses of the precompensated pattern is performed 
15 such that it assumes (almost) exclusively positive valu- 
es. A smooth solution is furthermore obtained since 
strong oscillations in the smearing function have no 
physical basis but are caused solely by mathematical 
instability of the calculations. 
20 In a preferred embodiment of the invention the 

method comprises the steps of : 

a) estimating a regularization parameter; 

b) determining a precompensated pattern with 
all pattern points of the desired pattern with the excep- 

25 tion of a determined pattern point; 

c) smearing the precompensated pattern again 
with the smearing function in order to predict the dose 
of the determined pattern point; 

d) repeating steps b and c for each pattern 

3 0 point ; 

e) repeating steps a to d with adapted regula- 
rization parameter until a final value of a regularizati- 
on parameter is obtained; 

f ) determining the precompensated pattern with 
35 the final value of the regularizat ion parameter. 

According to a further embodiment of the inven- 
tion step b) comprises the following iterative definiti- 
on : 
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with d (0) = 0 and r i0) = a 

in which a is a vector with the doses of the desired 
pattern as elements, d is a vector with the exposure 
5 doses of the precompensated pattern, K is the smearing 
function in matrix form, K" is the Hermitian conjugate of 
the smearing function K, B is an operator and X a 
regularization parameter. 

According to a further embodiment of the 
10 invention operator B is defined as follows: 

B(D) (-f^)ln(-^M 

i "cot ^toc 

in which the summation takes place over all pattern 
points, d. is the i ch element of the vector d, and d cot 
15 represents the summation over all elements of the vector 
d. 

According to a further embodiment of the 
invention the final value of the regularization parameter 
in the above mentioned step e) is the regularization 
20 parameter wherein 

in which N is the total number of pattern points, a is a 
. vector with the doses of the desired pattern as elements, 
d is a vector with the exposure doses of the 
25 precompensated pattern and K the smearing function in 
matrix form. 

According to a further embodiment of the 
invention the final value of the regularization parameter 
in step e) is the minimal regularization parameter 
3 0 wherein 
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±Y; <** - [Kd k (X) ] k ) 2 w kk ( X ) 

in which N is the total number of pattern points, a is a 
vector with the doses of the desired pattern as elements, 
d is a vector with the exposure doses of the 
precompensated pattern, K is the smearing function in 
5 matrix form and w kk is defined as: 

w kk a) = [ ]2 

with a kk the elements of the matrix A = K (K T K+A.L (D) t L (D) ) " 
l K T and L the discrete Laplace transformation. 

According to a further embodiment of the 
invention after step e) the step is performed of training 
10 a neural network using one or more desired first patterns 
and the associated precompensated patterns. 

According to a further embodiment of the 
invention the precompensated pattern associated with a 
second desired pattern can be determined with the trained 
15 neural network, wherein in a further embodiment the first 
desired pattern is a relatively simple training pattern 
and the second desired pattern is the partial pattern of 
an integrated circuit, and wherein in a further 
embodiment two or more partial patterns can be combined 
20 into a composite pattern of the integrated circuit. 

By determining the associated precompensated 
pattern of exposure doses for a known desired pattern, 
which is preferably simple, and then establishing the 
relation between the weighting factors of a neural 
25 network, is ensured that for a second desired pattern, 
which may be complicated, obtaining the relation between 
this pattern and the associated exposure doses is 
determined in very efficient and rapid manner by the 
neural network. The first pattern is aeneraily a 
30 relatively simple training pattern, while the second 
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pattern is for instance the pattern of a very complicated 
integrated circuit . 

In a preferred embodiment of the invention the 
above stated neural network is implemented in hardware, 
5 whereby determining of the relation between a pattern and 
the exposure dose associated therewith is performed in 
more rapid manner, for instance within 6 0 ns for each 
pattern point and within 10 minutes for a pattern of 10 10 . 

According to a preferred embodiment of the 
10 invention the smearing function is made up of at least 
two Gaussian functions, to which an exponential function 
is optionally added. Parameters of the Gaussian functions 
and optionally the exponential function can be determined 
by means of statistical simulation of the system of 
15 electron beam transmitting equipment and the relevant 
coating and substrate of the integrated circuit for 
manufacture . 

In another embodiment of the invention 
parameters are determined by performing measurements on 

2 0 the system of electron beam transmitting equipment and 

the relevant coating with substrate . 

The present invention also relates to a device 
for determining the exposure dose of an electron beam 
required per pattern position to obtain a desired pattern 
25 in a coating on a substrate, comprising an electronic 
circuit for implementing said neural network with 
weighting factors determined in the above stated manner. 

The invention will be elucidated hereinbelow 
with reference to a preferred embodiment thereof, wherein 

3 0 use will be made of the annexed drawings, in which: 

- figure 1 shows a schematic overview of a 
preferred embodiment of a device according to the 
invention; 

- figures 2a-2c show a schematic overview of 
35 the determination of a precompensated pattern of 3x3 

pattern points; 

- figure 3 shows a desired training pattern of 
256x256 pattern points; 



WO 99/66530 



PCT/BE99/00076 



- figure 4 shows the training pattern of figure 
3 after smearing; 

- figure 5 shows a graph in which for the 

training pattern of figure 3 the prediction error is 

5 plotted as a function of the chosen regularization 
parameter ,- 

- figure 6 shows the training pattern of figure 
3 after precompensation; 

- figure 7 shows the precompensated pattern of 
10 figure 6 after smearing; and 

- figure 8 is a schematic representation of a 
neural network for determining precompensated patterns. 

In an arrangement of equipment for transmitting 
an electron beam and a substrate 1 with coating 2 for 
15 processing, a beam of electrons 3 is directed at a 

position or pattern point of a coating 2 on a substrate 
1. The interaction of the incident electron beam 3 with 
the coating or resist film 2 and the underlayer or 
substrate 1 results in a scattering of the electrons in 
20 coating 2 which causes smearing or proximity effects. 

When for instance a primary electron penetrates into the 
coating, a part of its energy is transferred to electrons 
of the atoms of the coating, which causes ionization or 
excitation thereof. A collision between electrons with a 
25 large transfer of energy generates a secondary electron 
which generally has a direction of movement perpendicular 
to that of a primary electron. 

Smearing effects in electron beam lithography 
relate more generally to the process whereby the 
3 0 resolution of the exposed pattern is reduced by primary 
electron scattering (forward scattering) and secondary 
electron excitation (backward scattering) in the coating 
and the substrate of an integrated circuit for 
manufacture. Sharp features such as angles in the desired 
35 pattern are rounded, line thicknesses and interspaces are 
modified and in particular extreme cases some features 
even disappear completely or are merged in incorrect 
manner with adjacent features. 
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The smearing effects or proximity effects can 
be described by a smearing function, which shows the 
relation between on the one hand the exposure doses of a 
determined pattern point of a pattern for manufacturing 
5 in the coating and on the other the doses actually 
absorbed by this pattern point and adjacent pattern 
points. The effect of the smearing is thus established in 
the smearing function. 

Assuming that exposure and smearing are 

10 linearly and spatially invariant and that for a numeric - 
solution a discrete representation is preferred, the 
above can be expressed in matrix form as follows: A = KD, 
in which A is a column vector of which each element a. is 
the total energy dose which is actually absorbed in the 

15 associated pattern point, K is a smearing matrix of which 
each mn th element is the portion of the energy dose which 
is absorbed in pattern point m from a unit -exposure dose 
supplied to pattern point n, and D is a column vector 
made up of elements d ± which represent the exposure doses 

2 0 generated per pattern point by the electron beam 

equipment. Since the smearing effect is unavoidable, it 
is best to adapt the exposure doses d. of the different 
pattern points such that the dose a i actually absorbed in 
a pattern point is such that the desired pattern is still 
25 obtained. 

This so-called precompensation of the exposure 
dose of the electron beam can be performed according to 
the prior art by determining the inverse of the smearing 
matrix K. There are many ways of generally inverting a 

3 0 matrix. However, these methods often take no account of 

physical limitations, such as in this case for instance 
those of the electron beam transmitters. No negative 
exposure doses for instance are thus possible. A further 
drawback of such inversion methods is that the inverted 
35 matrix has many oscillations. In addition, for inversion 
of the smearing matrix for a partial pattern of for 
instance 256x256 pattern points the inversion of a 
smearing matrix with dimensions of 65536x65536 has to be 
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calculated, which requires an enormous amount of 
calculating time. 

Figures 2a to 2c show a desired pattern (A) . 
The pattern is built up of 9 pattern points a. wherein i 
5 varies from 1 to 9 . This desired pattern must' be 

precompensated in order to be able to provide the desired 
pattern after exposure to the smearing electron beam, 
i.e. the values of d ± , with i varying from 1 to 9, have to 
be determined. 
10 -~ The precompensated pattern is first of all 

determined making use of the doses ^ with i from 2 to 9 
wherein pattern point 1 is not therefore taken into 
account (figure 2a). This precompensated pattern is 
determined on the basis of the following expression: 

15 with d (0) = o and r c °> = a 

wherein a is a vector with the doses of the desired 
pattern as elements, d is a vector with the exposure 
doses of the precompensated pattern , K is the smearing 
function in matrix form, K* is the Hermitian conjugate of 

2 0 smearing function K, B is an operator and X a 

regularization parameter. The value of the regularization 
parameter can be chosen at random, in this case for 
instance A = 0 . 

The operator B imposes a limitation and can be 
25 defined as follows: 



B(D) (-^-)in(-^i-) 

i "coc a coc 



in which the summation takes place over all pattern 
points, d. is the i ch element of the vector d and d 

' cot 

represents the summation over all elements of the vector 



The thus determined precompensated pattern is 
then smeared once again on the basis of the known 
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smearing function, whereby the predicted dose Kd of 
pattern point 1 is determined. 

The above procedure is then repeated 
successively (figures 2b and 2c) for the second to ninth 
5 pattern point (i=2,..,9), wherein all pattern points with 
the exception of one pattern point are used each time. 

On the basis of the above results, the least 
squares prediction error over all pattern points is 
determined, which will be further explained later. 

10 The above procedure is subsequently repeated 

with different values for the regularizat ion parameter X, 
for instance A 2 = 0.001, A= 0.002 etc. The regularization 
parameter is eventually chosen wherein the. least squares 
prediction error over all pattern points is minimal. This 

15 regularization parameter is then chosen as the optimal 
regularization parameter A . The final precompensated 
pattern is then determined on the basis of this optimal 
regularization parameter A opt . 

For this purpose the minimum is determined of 

20 the expression: 

-if; (a k - [Ka fc (X)] Jt ) a w iHt a) 



in which N is the total number of pattern points, a is a 
vector with the doses of the desired pattern as elements, 
25 d is a vector with the exposure doses of the 

precompensated pattern, K is the smearing function in 
matrix form and w^ is defined as : 

WkkU ) = e 1 ~^* <X) ]» 

with a kk the elements of the matrix A = K ( K T K+ XL (D) T L (D) ) " * 
1 K T and L the Laplace operator. 
3 0 The smearing function resulting from forward 

scattering and backward scattering of the electrons of 
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the electron beam can be determined in different ways . it 
can be determined on the basis of measurements of the 
impulse" response of the equipment for transmitting the 
electron beam on a test obiect Thp • 
q a i=„ k ^ • J Tne shearing function can 

also be determxned using diverse Monte Carlo techniques, 
in the fxrst method of determination all physical aspects 
of the equipment used are taken into account, m the 
latter mentioned method of determination only a model of 
the reality is used, although the determination is 
10 however easier to perform without requiring extensive " 
measurements. 

Gaussian functions are preferably used as 
approximation for the smearing functions determined in 

15 abOVS deSCribSd methods " The smearing function 

15 ,n thxs case ..fitted" for instance with a scattering 

and T t °! 3 dOUblS GaUSSlan fUnCt±0n (f ° r bot * ^rwrd 
and backward scattering properties of the electrons) a 

trxple Gaussian function or a double Gaussian function 
wxth a decreasing exponential function. The choice of the 
20 scattering fit model depends on the dimensions of the 
components to be distinguished in the test object 
(resolution) . At dimensions smaller than 100 nm the 
choxce hereof becomes critical : at such small dimensions 
the trxple Gaussian functions or double Gaussian 
5 functions with decreasing exponential function are 
recommended. A smearing function with double Gaussian 
function can be described with 3 parameters, while the 
other two stated scattering fit models can be described 
wxth 4 parameters, which implies a great reduction in the 
0 quantity of data for processing. 

Figure 3 shows a desired pattern of 256x256 
pattern points. Smearing with a smearing function in the 
form of a double Gaussian function with «- 50 nm, S - 3 45 
and »» !.3 6 produces the smeared patt - ern Qf f . gure ^ ^ 

xs clearly visible that much detail in the pattern has 
been lost, which means a limitation in the resolution to 
be obtained of the pattern for manufacture. Application 
of the method according to the invention produces an 
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optimal regularization parameter of X opc = 0.07042, which is 
shown in' figure 5, in which the error in the pattern is 
minimal at this value of X. The precompensated pattern 
calculated with this value of X and the associated 

5 smeared pattern are shown respectively in figures 6 and 
7. Comparison of the results of figure 7 with those of 
figure 3 shows that the precompensation of the pattern 
with a desired pattern produces a smeared pattern with a 
greatly improved resolution. Components for 

0 distinguishing with dimensions of less than 100 nm, for - 
instance in an integrated circuit, can hereby be 
realized. A comparison of the results of the method 
described herein with those of other correction methods 
is shown in table 1. The degree of error of the 

5 correction methods is defined here as the summation of 
the difference between the calculated exposure doses and 
the ideal precompensated exposure doses divided by the 
number of pattern points . 



Correction method 


degree of error in % 


Uncorrected 


10 .2^ % 


Truncating 


10.2 % 


Shifting and scaling 


12.2 % 


. Present method 


4.9 % 



5 

From the above can be seen that the present method of 
determining a precompensated pattern produces by far the 
smallest degree of error compared to the other usual 
methods . 



0 The precompensated pattern and the desired 

pattern are subsequently used as training set or training 
patterns for a neural network. A part of such a network 
is shown schematically in figure 8 en is represented by 
the expression 
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i.e. the dose a. i is expressed in a set of 9 basic 
functions h ijt in this case radial functions. 

After training of the neural network a 
5 precompensated pattern can be determined for another 
random desired pattern in very rapid manner. A random 
pattern can for instance be a pattern of 512 by 512 
pattern points forming a partial pattern of an integrated 
circuit. Various partial patterns can then be combined 
10 (clustered) to form one pattern which comprises the whole 
integrated circuit or at least a part thereof. 

The above described neural network can be 
implemented in hardware, and preferably in analog 
hardware since the calculating speed of neural networks 
15 implemented in this manner is very great. The calculating 
time for precompensation of a pattern thus amounts to 
less than 60 ns per pattern point. Precompensation of a 
pattern of an integrated circuit of about 10 10 pattern 
points requires in this case only fcbout 10 minutes on 
20 present personal computers. 

The invention is further described in the non- 
prepublished doctoral thesis with the title "Proximity 
effects correction in electron beam nanolithography" , the 
entire content of which should be deemed as interpolated 
25 herein. 



